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ABST21ACT 


Under  certain  meteorological  conditions  ships  may  generate  "tracks"  in  overlying  marine 
stratus  or  stratocumulus  clouds.  These  tracks  are  often  not  particularly  evident  to  satellite  sensors 
at  visible  wavelengths,  but  show  up  much  more  clearly  at  near  infrared  wavelengths,  such  as 
chaimel  3  (3.7  pm)  of  the  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  on  NOAA 
polar-orbiting  satellites.  This  phenomenon  represents  a  potential  naval  intelligence  technique,  not 
only  as  a  means  to  monitor  and  perhaps  classify  ships  of  other  nations  but  also  as  a  measure  of 
envirorunental  vulnerability  of  U.S.  ships.  A  simple  theoretical  model  was  used  as  part  of  a  ship 
track  study  to  investigate  effects  of  changes  in  microphysics  and  geometry  on  cloud  reflectance 
at  certain  solar  wavelengths.  The  model  assumed  a  plane  parallel  homogeneous  atmosphere  of 
infinite  extent  and  finite  thickness,  with  Mie  scattering  and  the  delta-Eddington  approximation 
for  radiative  transfer.  The  results  were  consistent  with  the  theory  that  ship  exhaust  can  increase 
the  number  density  of  smaller  cloud  droplets,  which  would  lead  to  an  increase  in  reflectance  of 
incident  solar  radiation,  especially  at  wavelengths  near  3.7  pm. 
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A  SIMPLE  CLOUD  REFLECTANCE  MODEL  FOR  SHIP  TRACKS  IN  CLOUDS 


1.  INTRODUCTION 

Clouds  are  a  ubiquitous  phenomenon,  typically  covering  about  fifty  percent  of  the  globe. 
They  play  a  crucial  role  in  determining  the  radiation  budget  of  the  earth  and  hence  have  a  major 
effect  on  the  dynamics  of  the  atmosphere  and  the  earth’s  climate.  It  has  also  been  established 
that  clouds  are  sensitive  to  anthropogenic  effects,  such  as  pollution,  thus  providing  a  link  between 
man’s  activities  and  his  environment.  One  intriguing  naval  phenomenon  is  the  generation  of 
"ship  tracks"  in  clouds  (Siquig  et  al.,  1991).  An  example  is  given  in  Figure  1.  which  is  a  channel 
3  (3.7  pm)  image  from  the  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  on  NOAA 
polar  orbiter  satellites  showing  a  number  of  ship  tracks.  In  AVHRR  channel  1  (.63  pm),  the 
tracks  are  much  less  evident.  One  possible  explanation  for  ship  tracks  is  that  ship  exhaust  can 
increase  the  number  of  cloud  condensation  nuclei  in  overlying  marine  stratus  or  stratocumulus 
clouds,  resulting  in  more  but  smaller  cloud  droplets  which  increase  the  cloud  reflectance  at  near 
infrared  (NIR)  wavelengths  relatively  more  than  at  visible  (VIS)  wavelengths  (Coakley  et  al., 
1987).  Here  a  simple  cloud  reflectance  model  is  used  to  explore  this  suggestion. 

We  start  by  observing  that  ship  tracks  usually  occur  in  marine  stratus  or  stratocumulus 
(st/sc)  cloud  decks  which  usually  are  some  fraction  of  a  kilometer  thick.  The  observed  tracks 
themselves  seem  to  range  in  width  from  a  few  to  about  twenty  kilometers.  As  a  first 
approximation,  then,  a  ship  track  can  be  considered  as  an  infinite  plane  parallel  homogeneous 
st/sc  cloud  layer  with  specified  microphysics  and  thickness.  This  allows  us  to  follow  the 
approach  of  Mineart  (1988),  hereafter  "Mineart,"  who  demonstrated  for  marine  stratocumulus 
clouds  the  sensitivity  of  reflectance  to  cloud  microphysics  and  geometry.  We  have  extended 
Mineart’s  analysis  to  include  more  data  points  and  additional  wavelengths,  resulting  in  an 
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increase  in  the  data  base  of  theoretical  results  by  a  factor  of  ten.  Results  were  in  basic  agreement 
with  the  earlier  work.  The  exhaust  explanation  for  ship  tracks  is  consistent  with  the  model 
results,  which  show  that  a  shift  to  smaller  cloud  droplet  size  has  a  greater  impact  proportionally 
on  cloud  reflectance  at  NIR  than  at  VIS  wavelengths. 

Section  2  is  a  brief  review  of  relevant  Mie  scattering  theory,  including  scattering  from  a 
single  homogeneous  sphere  and  scattering  from  a  polydispersion  of  homogeneous  spheres  with 
a  specified  size  distribution.  Section  3  reviews  the  delta-Eddington  approximation  to  the  radiative 
transfer  equation  for  a  plane-parallel  homogeneous  cloud  layer  of  infinite  extent.  The  choice  of 
parameter  values  is  described  in  Section  4.  Representative  results  for  reflectance  as  a  function 
of  several  variables  (wavelength,  size  distribution,  mode  radius,  equivalent  liquid  water  content, 
cloud  thickness,  and  solar  zenith  angle)  are  presented  in  Section  5.  Conclusions  and  suggestions 
for  fuiure  work  are  given  in  Section  6. 

2.  RELEVANT  MIE  SCATTERING  THEORY 

2. 1  Overview 

Attenuation  of  radiation  in  the  atmosphere  occurs  principally  through  absorption  and 
scattering  by  atmospheric  gases,  cloud  droplets  and  paniculate  maner.  By  confining  our  attention 
to  atmospheric  windows,  the  problem  of  atmospheric  gases  can  be  largely  avoided,  and  for  the 
wavelengths  chosen,  aerosols  are  not  as  imponant  as  cloud  droplets  because  of  their  smaller 
scattering  efficiencies.  Thus,  this  simple  model  is  reduced  to  a  problem  in  cloud  attenuation. 
To  sinqjlify  the  problem  even  further,  we  assume  a  plane  parallel  homogeneous  cloud  of  infinite 
extent.  The  problem  is  then  separated  into  several  parts.  In  Section  2.2,  scattering  of 
electromagnetic  (EM)  unpolarized  plane  waves  off  a  single  homogeneous  sphere  is  reviewed. 
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The  results  are  extended  in  Section  2.3  to  a  polydispersion  of  homogeneous  spheres  with  a 
specified  droplet  size  distribution. 

2.2  Single  Particle  Mie  Scattering 

Maxwell’s  equations  can  be  used  to  generate  a  vector  wave  equation  for  electromagnetic 
waves  in  a  given  medium.  For  a  homogeneous  sphere  with  complex  index  of  refraction  m, 

m  =  v-ik 

where  the  real  and  imaginary  parts  have  a  wavelength  dependence,  the  incident  and  scattered 
wave  solutions  can  be  matched  at  the  sphere  boundary  to  arrive  at  the  Mie  scattering  parameters 
for  the  far  field.  Stokes  vectors  (related  to  intensities)  for  the  incident  and  scattered  fields  are 
related  by  a  transformation  matrix  from  which  the  scattering  phase  matrix  is  derived.  A 
prominent  feature  of  Mie  scattering  is  its  strong  forward  bias,  a  feature  which  will  be 
fundamental  for  explaining  many  of  the  computational  results. 

It  is  convenient  to  introduce  a  dimensionless  size  parameter  x,  which  is  the  ratio  of  the 
sphere’s  circumference  to  the  radiation  wavelength, 

x  =  2nrlX  ^2) 


where 


and 


r  =  sphere  radius 


A.  =  wavelength. 


For  the  two  extremes  r  «  X  (Rayleigh  scattering)  and  r  »  X  (geometric  optics),  simple 
scattering  formulae  apply.  For  the  intermediate  regime,  r  of  the  order  of  X,  the  correct  solution 
(Mie)  is  quite  complicated  and  is  generally  calculated  only  through  computer  programs.  Since 
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cloud  droplet  sizes  are  comparable  to  the  NIR  wavelengths  of  interest  (micrometer  range),  it  was 
necessary  to  use  such  a  program  for  this  study  (see,  e.g.,  Wiscombe,  1980). 

The  three  basic  quantities  that  describe  single  panicle  scattering  are 
Q,  (m,r)  =  scattering  efficiency  factor 
Q,  {mjc)  =  extinction  efficiency  factor 
P  (m,r,\j/)  =  phase  function  for  radius  r, 

where  y  is  the  scattering  angle  relative  to  the  incident  direction.  Basically,  the  efficiency  factors 
are  a  measure  of  process  cross-section  of  the  particle  relative  to  its  geometric  cross-section,  while 
the  phase  function  is  a  measure  of  the  directionality  of  the  scattering  process.  Figure  2  shows 
the  complicated  nature  of  Q„  depending  on  the  value  of  k. 


5  10  50  100 

X  -Ztr  r/\ 


Figure  2.  Mie  scattering  efficiency  factor  (Q,):  Illustrated  as  a  function  of  the  size  parameter 
X  =  litr/X.  The  real  refractive  index  is  1.33,  with  results  shown  for  four  values  of 
the  imaginary  pan  of  the  refractive  index  (after  Hansen  and  Travis,  1974). 
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2.3  Polydispersion  Mie  Scattering 

The  Mie  scattering  analysis  must  be  extended  to  a  collection  of  homogeneous  spheres  in 
order  to  be  applicable  to  a  cloud  model.  Consider  a  polydispersion,  i.e.,  a  colleaion  of  scattering 
spheres  of  uniform  physical  properties,  differing  only  in  their  size  according  to  a  known 
continuous  size  distribution.  Furthermore,  assume  the  particles  are  separated  enough  so  that 
scattering  off  each  particle  is  independent  of  scattering  off  the  other  particles.  Then  the 
intensities  of  scattering  by  various  particles  can  be  added  without  considering  the  phases  of  the 
scattered  waves  ("independent  scattering"). 

By  specifying  a  droplet  size  distribution  n(r)  it  becomes  possible  to  generalize  the  above 
parameters  to  those  appropriate  for  a  polydispersion: 
volume  scattering  coefficient 

o,(A)  =yj‘7tr^n(r)<7,(m,r)d'r  (3) 


volume  extinction  coefficient 

o,(^)  =  f’  itr^n(r)Q^im,r)dr 

*  O 


(4) 


scattering  phase  function 

=  —77-  r  nr^n(r)Q  (m,r)P(m,r,ifr)dr .  (5) 
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In  addition  two  other  related  parameters  can  be  defined,  the  single  scattering  albedo 


o.(A) 


(6) 


which  is  related  to  the  probability  that  scattering  will  occur  instead  of  absorption,  and  the 
asymmetry  factor  g 


g  -  <cos\|r>  =  1/2  P,(m,4r)  cos\|k/cos<1; 


(7) 


which  is  related  to  the  first  moment  of  the  phase  function. 

The  actual  details  of  the  Mie  derivation  of  Q,  (w,r),  (m,r)  and  P(m,r,y)  are  available 
in  standard  references,  e.g.,  Kerker  (1969)  and  Liou  (1980);  suffice  it  to  say  that  the  required 
parameters  are  x,  m,  and  X.  With  the  specification  of  n(r)  instead  of  x,  then  the  volume 
scattering  coefficient  the  single  scattering  albedo  (o^,  and  the  asymmmetry  factor  g  can  be 
used  to  characterize  scattering  off  a  small  volume  element  of  homogeneous  spheres  with  size 
distribution  n(r).  With  suitable  simplifying  approximations,  these  three  parameters  are  sufficient 
to  describe  the  transfer  of  radiation  within  a  cloud. 

Before  discussing  in  Section  3  the  radiative  transfer  treatment,  it  might  be  instructive  to 
point  out  why  these  three  parameters  are  so  useful.  The  volume  scatteiing  coefficient  leads 
directly  to  the  scattering  optical  depth  S^(X,z)  for  a  physical  distance  z, 


which  for  a  homogeneous  medium  is  equivalent  to 

5,(A,z)  =  a/A)z 


(8) 


(9) 
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where  5,  is  an  optical  measure  of  attenuation  from  scattering.  The  same  physical  distance  z  can 
have  different  values  for  depending  on  a,.  This  fact  will  be  useful  later  in  interpreting  the 
model  results.  The  single  scattering  albedo  is  a  measure  of  scattering  versus  absorption  in  the 
interaction  of  incident  radiation  with  a  small  volume  element.  For  example,  to„  =  1  means  pure 
scattering  with  no  absorption.  Finally,  the  asymmetry  factor  is  an  integrated  measure  of  the 
scattering  directionality,  where  g  =  1  means  complete  forward  scattering,  g  =  0  means  isotropic 
scattering,  and  g  =  -1  means  complete  backward  scattering.  Note  that  if  an  effective  radius  r^ 
is  defined  for  the  distribution  such  that 


r 


L 

/: 


r^n{f)dr 


rhi{r)dr 


(10) 


and  if  we  define  x, 

x^  =  2nrJX  (11) 

then  g  approaches  0  as  x,  becomes  small  (Rayleigh  regime)  and  g  approaches  approximately  .87 
as  X,  becomes  large  (geometric  optics). 

Thus,  under  certain  conditions,  these  three  parameters  are  sufficient  to  describe  the 
integrated  effects  of  all  the  radiation  scattered  from  an  extended  cloud  layer  of  finite  thickness. 


3.  RADIATIVE  TRANSFER  TREATMENT 

3.1  Multqjle  Scattering  and  Radiative  Transfer 

Most  water  clouds  are  optically  thick,  that  is  very  few  incident  photons  will  escape 
scattering.  Since  the  probability  for  multiple  scattering  is  very  high  in  a  cloud,  then  the  details 
of  the  phase  function  for  single  particle  scattering  are  smoothed  out  in  the  multiple  scattering 
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process,  and  only  gross  features  of  the  phase  function  need  to  be  specified,  such  as  the 
asymmetry  factor.  Note,  though,  that  if  absorption  becomes  important  through  the  wavelength 
dependence  of  the  index  of  refraction,  then  the  resultant  scattered  radiation  retains  more  of  its 
directional  cheu'acteristics.  This  is  because  any  scattered  radiation  cannot  have  entered  too  deeply 
into  the  cloud  and  still  emerged,  so  it  carmot  have  undergone  too  many  scatterings.  This 
connection  between  directionality  of  scattered  radiation  and  the  absorption  properties  of  the  cloud 
will  be  useful  in  explaining  some  of  the  model  results. 

Recall  that  a*,  (Og,  and  P,  respectively  determine  1)  how  much  energy  is  lost  from  the 
incident  beam  in  traversing  a  small  volume  element,  2)  the  probability  that  this  energy  will  be 
lost  by  scattering  or  absorption,  and  3)  the  probability  that  scattered  energy  will  be  sent  into  a 
particular  direction.  For  a  cloud  composed  of  similar  small  volume  elements,  its  effect  upon 
incident  radiation  will  also  depend  upon  the  geometry  of  the  cloud  and  its  orientation  with 
respect  to  the  incident  beam.  For  hemispheric  reflectance,  the  solar  zenith  angle  will  provide  the 
cloud-beam  orientation.  The  geometry  problem  can  be  simplified  by  choosing  a  homogeneous 
plane  parallel  layer  of  infinite  extent.  These  assumptions  reduce  the  complexity  of  the  radiative 
transfer  equation  enormously  and  allow  us  to  use  the  delta-Eddington  approximation,  described 
below.  So  our  additional  parameters  are  cloud  thickness  z  and  solar  zenith  angle  (sza)  O^. 

3.2  Delta-Eddington  Approximation 

Most  solutions  to  specific  problems  in  radiative  transfer  require  complex  computer 
implementation.  For  our  simple  model  a  major  simplification  can  be  made,  namely  that  the 
radiance  can  be  described  as  a  simple  two-term  expansion  (Eddington  approximation  -  see  Shettle 
and  Weinman.  1970).  This  results  in  straight-forward  expressions  for  the  reflectance  off  an 
infinite  homogeneous  plane-parallel  cloud  layer.  The  only  input  required  is  solar  zenith  angle, 
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asymmetry  factor,  scattering  optical  depth  and  single  scattering  albedo.  This  approximation  can 
be  improved  by  exploiting  the  significant  forward  peak  in  the  scattering  phase  function.  If  the 
forward  peak  is  approximated  by  a  Dirac -delta  function  and  the  phase  function  by  a  two  term 
expansion,  then  the  resultant  delta-Eddington  approximation  is  equivalent  to  the  Eddington 
approximation  with  transformed  co,,,  g,  and  5,  (Joseph  et  al.,  1976).  Thus  the  simple  Eddington 
expressions  for  reflectance  can  be  used  with  slight  modification,  vastly  reducing  the 
computational  burden. 

3.3  Additional  Comments 

We  have  seen  that  specifying  the  droplet  composition  m,  size  distribution  n(r),  and 
wavelength  X,  allow  us  to  determine  the  scattering  optical  depth  5,.  single  scattering  albedo 
and  asymmetry  factor  g.  These,  in  turn,  allow  us  to  calculate  hemispheric  reflectance,  given 
cloud  thickness  and  solar  zenith  angle.  We  have  assumed  perfect  transmission  above  the  infinite 
cloud  layer  because  of  the  window  regions  in  the  EM  spectmm  we  are  using.  The  attenuation 
coefficient  of  aerosols  above  the  cloud  is  small  compared  to  that  for  the  cloud  droplets,  and  also 
the  water  vapor  content  above  the  cloud  dies  off  rapidly. 

4.  CHOICE  OF  PARAMETER  VALUES 

4.1  Approach 

Our  interest  is  in  the  effect  of  cloud  microphysics  (droplet  size  distribution  and  equivalent 
liquid  water  content  (ELWC))  and  geometry  (cloud  thickness)  on  measured  reflectance  values 
in  atmospheric  transmittance  windows.  To  facilitate  comparison  with  Mineart’s  work  and  extend 
its  analysis,  we  have  used  fiis  choice  of  mode  radii,  distributions,  wavelengths,  ELWC,  cloud 
thicknesses,  and  solar  zenith  angles  while  adding  extra  values. 


10 


4.2  Parameter  Values 


The  following  wavelengths  were  chosen;  .63,  1.6,  2.3,  and  3.7  pm.  The  first  and  last  are 
channels  on  the  AVHRR  and  were  used  by  Mineart.  The  other  two  are  possible  channels  on 
future  satellites.  All  these  wavelengths  are  in  atmospheric  windows. 

A  modified  gamma  distribution  as  discussed  in  Deirmendjian  (1969)  was  used  to  model 
droplet  size  distribution. 


«(r)=Ar*c 


(12) 


where  A,  a,  y.  tc  positive  and  real  and  a  is  an  integer.  Note  that  n(0)  =  0  and  n(«>)  =  0  and 
that  n(r)  =  maximum  at  r  =  r^,  the  mode  radius.  The  distribution  n(r)  gives  the  number  density 
of  particles  within  a  small  radius  range.  The  total  number  density  of  panicles  is  given  by  N, 

N  =  C  nir)dr ,  U3) 

J  O 


from  which  it  is  clear  that  A  depends  on  N.  The  three  distributions  used  by  Mineart  were 


DISTRIBUTION 

a 
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D1 

2 

1.19 

D2 

5 

2.41 

D3 

5 

1.30 

The  values  chosen  for  r^  were  4,  6,  8,  10,  and  12  pm  to  bracket  expected  values  for  marine 
stratocumulus  clouds,  for  which  in  situ  droplet  data  were  available.  Mineart  used  the  values  4 
and  8  for  r^.  For  ELWC  the  following  values  were  used:  .2,  .4,  .6,  and  .8  gm/m^  Mineart  used 
the  values  .4  and  .8. 
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By  specifying  ELWC  it  was  necessary  to  calculate  the  corresponding  values  for  A.  This 
was  done  by  recognizing  that 


ELWC  =  5V  =E 


(14) 


where 


6  =  density  of  water  =10^  gm/m’ 
V  =  total  droplet  volume 


4/3  nr^n{r)dr . 


(15) 


It  is  straightforward  to  show  that 

g  *4 

Y(--)  " 

A  = - 1 - 

4/37:r;*"r(-55-^)p 

Y 

where  F  is  the  normal  gamma  function.  Tj'pical  curves  for  Dl,  D2,  and  D3  are  shown  in  Figure 
3,  taken  from  Mineart.  D3  may  be  considered  as  a  typical  droplet  size  distribution  for  marine 
stratocumulus  clouds,  D2  is  a  narrow  peaked  distribution,  and  Dl  is  a  broad  distribution.  These 
distributions  will  be  used  to  examine  the  effect  of  distribution  shape  on  reflectance  values. 

The  values  of  cloud  thickness  chosen  were  those  used  by  Mineart :  10,  25,  50,  100,  250, 
500,  and  750  m.  The  solar  zenith  angle  values  used  were  also  those  of  Mineart :  30, 45,  and  60 
degrees. 

Figure  4  is  a  schematic  of  the  scattering  process  being  modeled. 
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5. 


RESULTS 


5.1  Overview 

The  resulting  data  base  of  hemispheric  reflectance  as  a  function  of  wavelength, 
distribution,  mode  radius,  ELWC,  cloud  thickness,  and  solar  zenith  angle  consisted  of  5040 
entries,  an  increase  by  a  factor  of  ten  of  Mineart’s  published  results.  Representative  results  will 
be  presented  here. 

The  addition  of  two  more  wavelengths,  three  more  mode  radii,  and  two  more  values  for 
ELWC  provide  a  ntore  complete  basis  for  asserting  that  an  NIR  channel  could  provide  an  indica¬ 
tor  of  cloud  microphysical  conditions.  Furthermore,  wavelengths  at  1.6  and  2.3  pm  also  look 
promising  for  this  purpose  in  addition  to  3.7  pm  as  identified  by  Mineart.  The  shorter 
wavelengths  have  the  advantage  of  greater  energy  in  the  incident  solar  spectrum  and  higher 
reflectance  values,  thus  enhancing  satellite  sensor  performance.  Combining  the  results  of  several 
NIR  channels  may  provide  a  more  powerful  tool  than  a  single  channel  alone.  Consistent  with 
the  earlier  results,  VIS  wavelengths  (.63  pm)  are  more  suited  to  characterizing  cloud  geometry 
(thickness)  than  microphysics. 

5.2  Scattering  Optical  Depth  vs.  Cloud  Thickness 

Figure  5  shows  that  for  this  homogeneous  model  scattering  optical  depth  increases  linearly 
with  thickness  as  expected.  The  effect  of  increasing  mode  radius  is  to  decrease  optical  depth. 
This  may  be  explained  by  considering  the  effective  radius  for  the  distribution,  r,.  The  scattering 
optical  depth  has  the  dependence 


(17) 
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Figure  5.  Scattering  Optical  Depth  (5j  vs.  Qoud  Thickness  (z).  For  ELWC  =  .4  gm/m^  and 
the  D3  distribution.  The  four  families  of  curves  are  labeled  by  mode  radius  and  within  each 
family  the  curves  are  distinguished  by  wavelength. 


where  for  a  fixed  ELWC 


ELWC 

p4/37irj 


(18) 


Therefore, 


indicating  that  if  effective  radius  alone  is  increased,  then  scattering  optical  depth  will  decrease. 
Note  that  it  is  straightforward  to  show  that 


(20) 


so  that  r,  and  r^  are  related  by  a  constant  for  a  given  distribution,  so  the  above  result  applies  also 
for  the  mode  radius. 

Increasing  the  wavelength  may  either  increase  or  decrease  the  slope  of  the  optical  depth 
vs.  thickness  curve,  depending  on  the  mode  radius.  For  our  homogeneous  approximation,  optical 
depth  is  linearly  related  to  the  volume  scattering  coefficient  given  by  equation  3;  therefore,  this 
effea  is  probably  due  to  the  oscillatory  nature  of  the  Mie  scattering  efficiency  factor  as  a 
function  of  x  and  index  of  refraction  (which  has  a  wavelength  dependence). 

From  equation  16,  we  see  that  A  scales  as  ELWC  and  therefore  so  docs  a,  and  hence  also 

8.- 

5.3  Single  Scattering  Albedo  vs.  Wavelength 

Figure  6  shows  that  as  wavelength  increases  the  single  scattering  albedo  decreases.  This 
is  a  consequence  of  the  generally  greater  absorption  at  the  longer  wavelengths,  which  is  a  result 
of  the  larger  imaginary  part  of  the  index  of  refraction  for  water  at  the  longer  wavelengths.  For 
a  given  wavelength,  as  mode  radius  increases,  single  scattering  albedo  decreases  (if  the 
wavelength  is  long  enough).  Note  that  is  more  sensitive  to  changes  in  r^  at  the  longer 
wavelengths. 
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5.4  Asymmetry  Factor  vs.  Wavelength 

Figure  7  shows  that  g  does  not  have  a  simple  monotonic  behavior  with  respect  to 
wavelength,  dep)ending  on  mode  radius.  This  makes  it  harder  to  generalize  about  the 
computational  results  on  the  basis  of  the  behavior  of  g.  There  does  seem  to  be  a  tendency  for 
g  to  decrease  as  wavelength  is  increased,  i.e.,  as  absorption  is  increased.  For  a  given  wavelength, 
g  increases  as  mode  radius  increases,  indicating  increased  forward  scattering.  This  effect  is  more 
pronounced  for  the  longer  wavelengths. 

5.5  Reflectance  vs.  Mode  Radius 

Figures  8a  through  8d  show  clearly  that  1)  reflectance  decreases  as  mode  radius  increases; 
2)  reflectance  decreases  as  wavelength  increases;  3)  reflectance  increases  as  ELWC  increases, 
with  this  effect  being  less  at  the  longer  wavelengths;  and  4)  reflectance  increases  as  cloud 
thickness  increases,  with  this  effect  being  less  at  the  longer  wavelengths.  The  laner  two  points 
are  better  described  in  later  figures.  For  a  given  ELWC,  as  r^  increases,  the  number  of 
droplets/unit  volume  must  decrease.  Therefore,  for  a  given  ELWC  fewer  but  larger  particles  have 
less  reflectance  than  more  but  smaller  particles.  One  reason  is  that  the  larger  particles  tend  to 
have  more  forward  and  less  backward  scattering.  The  second  point  above  is  primarily  due  to  the 
increased  absorption  at  the  longer  wavelengths. 

5.6  Reflectance  vs.  ELWC 

Figures  9a  through  9d  show  that  1 )  in  general  reflectance  increases  as  ELWC  increases; 
2)  this  effect  is  less  pronounced  at  the  longer  wavelengths;  3)  at  any  wavelength  reflectance 
seems  to  saturate  beyond  a  certain  value  for  ELWC;  4)  reflectance  decreases  as  wavelength 
increases;  5)  reflectance  increases  as  cloud  thickness  increases,  up  to  a  point;  and  6)  reflectance 
decreases  as  mode  radius  increases.  The  first  point  simply  means  that  at  fixed  r,.  an  increase  in 
ELWC  means  an  increase  in  the  number  of  particles,  and  therefore  more  incident  radiation  is 
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Figure  8.  Reflectance  vs.  Mode  Radius.  For  all  cases  the  distribution  is  D3  and  =  45°. 
ELWC  and  cloud  thickness  vary  as  indicated.  The  curve  symbols  for  the  different  wavelength 
values  are  used  consistently  in  the  following  figures. 
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SLWC  (gm/m*)  ELWC  (gm/m*) 


Figure  9.  Reflectance  vs.  ELWC.  For  all  cases  the  distribution  is  D3  and  0„  =  45°.  Mode 
radius  and  cloud  thickness  vary  as  indicated. 
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reflected.  The  second  point  is  due  to  the  increased  importance  of  absorption  at  the  longer 
wavelengths.  The  third  and  fifth  points  are  best  discussed  in  terms  of  optical  depth.  When  the 
optical  depth  of  the  cloud  exceeds  a  certain  value,  it  may  as  well  be  infinite;  increasing  it  no 
longer  has  any  effect  on  reflectance.  Optical  depth  (recall  Figure  5)  can  depend  on  physical 
thickness,  wavelength  (through  the  effect  of  absorption  on  a.)  and  ELWC.  As  previously 
discussed  the  fourth  point  can  be  explained  by  the  increased  absorption  at  the  longer  wavelengths. 
The  last  point  was  explained  in  the  first  set  of  figures. 

5.7  Reflectance  vs.  Qoud  Thickness 

Figures  10a  through  lOd  show  that  1)  up  to  a  certain  point  reflectance  increases  as 
thickness  increases;  2)  reflectance  decreases  as  wavelength  increases;  3)  reflectance  decreases  as 
mode  radius  increases;  4)  up  to  a  point  reflectance  increases  as  ELWC  increases;  and  5)  the 
cloud  thickness  at  which  reflectance  saturates  decreases  as  wavelength  increases  or  as  mode 
radius  decreases.  The  first  and  fourth  points  have  been  explained  from  the  point  of  view  of 
optical  depth.  The  second  point  has  been  explained  by  greater  absorption  at  longer  wavelengths. 
The  third  point  is  due  to  fewer  scatterers  at  a  given  ELWC.  The  first  part  of  the  fifth  point  is 
due  to  greater  absorption  at  longer  wavelengths,  allowing  a  thinner  cloud  to  reach  saturation  for 
reflectance.  The  second  part  of  the  fifth  point  is  because  more  scatterers/unit  volume  means  an 
increase  in  optical  depth  for  a  given  physical  thickness. 

5.8  Dependence  of  Reflectance  on  Distribution 

Figure  1 1  shows  that  reflectance  is  greatest  for  the  D2  distribution,  followed  by  D3  and 
then  Dl.  D2  is  a  narrow  peaked  distribution,  D1  has  a  broad  peak,  and  D3  is  intermediate. 
From  equation  20  one  can  calculate  the  ratio  of  effective  radius  to  mode  radius.  For  the  D2,  D3, 
and  Dl  distributions  this  ratio  is  1.17,  1.42,  and  2.13  respectively.  Thus,  the  dependence  of 
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Reflectance 


1.00 


reflectance  on  distribution,  assuming  the  same  mode  radius  and  ELWC,  essentially  is  because 
more  and  smaller  droplets  have  higher  reflectance. 

5.9  Dependence  of  Reflectance  on  Solar  Zenith  Angle 

Figures  12  through  14  show  that  reflectance  increases  as  solar  zenith  angle  increases  over 
the  range  considered,  and  that  this  effect  seems  to  be  proportionally  more  pronounced  for  the 
longer  wavelengths.  The  first  point  can  be  explained  by  the  strong  forward  scattering 
characteristic  of  this  Mie  size  range.  Increasing  sza  allows  more  scattering  into  the  upper 
hemisphere.  The  second  point  is  mainly  because  for  the  longer  wavelengths  the  scattered 
radiation  is  more  anisotropic  since  it  sees  a  thinner  layer  of  the  cloud  due  to  the  greater 
absorption  at  longer  wavelengths. 

6.  CONCLUSIONS,  FUTURE  WORK 

6. 1  Conclusions 

The  calculated  reflectance  values  depend  on  a  number  of  different  variables,  not  all 
independent.  Some  of  the  dependencies  are  simple,  while  others  are  more  con^licated.  For 
example,  the  complex  nature  of  the  index  of  refraction  as  a  function  of  wavelength  makes 
interpretation  of  the  reflectance  values  in  terms  of  absorption  properties  of  water  subject  to 
qualifications.  Nonetheless,  if  only  one  variable  at  a  time  is  varied,  its  effect  on  reflectance  can 
be  determined,  at  least  over  the  range  of  values  considered  here.  From  the  calculations  one  can 
then  derive  the  following  tendencies: 
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SUMMARY  OF  VARIATIONS 


VARIABLE 

CHANGE 

REFLECTANCE 

increase 

decrease 

wavelength 

increase 

decrease 

thickness 

increase 

increase 

ELWC 

increase 

increase 

sza 

increase 

increase 

width(r,) 

increase 

decrease 

As  implied  above,  some  of  these  changes  do  not  occur  without  limit.  A  few  remarks  are 
in  order  to  clarify  these  dependencies. 

1.  When  conditions  combine  to  give  the  cloud  a  sufficiently  large  optical  depth 
(essentially  infinite),  then  any  change  that  acts  to  increase  optical  depth  has  little  effect  on 
reflectance.  For  example,  increasing  thickness  or  ELWC  beyond  certain  limits  no  longer  affects 
reflectance. 

2.  For  a  given  wavelength  and  solar  zenith  angle,  the  critical  cloud  factors 
determining  reflectance  are  single  scattering  albedo,  asymmetry  factor,  and  optical  depth.  As 
wavelength  varies,  the  values  of  these  factors  wiU  change  and  as  a  consequence  their  relative 
importance  will  change.  Thus,  an  effect  which  may  be  dominant  for  one  set  of  values  may  not 
be  critical  for  another  set. 

3.  The  importance  of  the  asymmetry  factor  and  the  single  scattering  albedo 
ultimately  depends  on  the  directional  properties  of  Mie  scattering,  which  are  determined  from  the 
relative  size  of  the  scatterers  compared  to  the  wavelength  of  the  radiation,  and  on  the  wavelength 
dependence  of  the  complex  index  of  refraction  of  water. 

4.  If  ELWC  is  kept  constant,  then  decreasing  the  effective  radius  of  the  droplets 
results  in  more  and  smaller  droplets  which  increase  reflectance.  This  is  relevant  to  the  ship  track 
phenomenon,  as  discussed  below. 
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5.  Switching  to  a  longer  wavelength  reduces  the  reflectance  in  general  because  of 
increased  absorption.  Note  that  this  is  based  on  the  results  for  four  wavelengths.  Because  of  the 
undulatory  nature  of  the  imaginary  part  of  the  index  of  refraction,  this  observation  may  not  apply 
at  all  wavelengths. 

6.  Subject  to  the  qualification  discussed  in  point  1,  an  increase  in  ELWC  will 
increase  reflectance  until  optical  depth  saturates.  A  similar  conclusion  holds  for  cloud  thickness. 
Larger  particles  require  a  larger  value  of  ELWC  for  optical  depth  to  saturate,  i.e.,  for  reflectance 
to  become  insensitive  to  further  increase  in  ELWC.  Similarly,  for  larger  particles  a  larger 
thickness  is  required  for  reflectance  to  saturate.  The  reasons  for  this  were  discussed  in  Section 
5.2. 

7.  Variations  in  solar  zenith  angle  have  a  larger  effect  on  reflectance  for  longer 
wavelengths. 

8.  Effects  of  distribution  can  be  explained  mainly  in  terms  of  differing  effective 
radius,  as  mentioned  in  point  4. 

From  the  above  results  two  general  recommendations  can  be  made  regarding  wavelength 
selection  for  determining  cloud  conditions: 

1.  The  shorter  wavelengths  (e.g.,  .63  pm)  are  more  useful  for  determining  cloud 
thickness  because  the  longer  wavelengths  see  a  thinner  layer  of  the  cloud  due  to  absorption 
limitations. 

2.  On  the  other  hand,  the  longer  wavelengths  are  more  sensitive  to  the  cloud 
microphysics,  such  as  droplet  size  and  ELWC.  Since  for  NIR  solar  radiation,  incident  energy 
drops  off  with  increasing  wavelength  and  reflectance  decreases  at  the  longer  wavelengths,  then 
the  shortest  wavelength  that  remains  sensitive  to  cloud  microphysics  should  be  used.  This  has 
tfie  advantage  of  increasing  receiver  sensitivity.  Mineart  had  previously  shown  that  3.7  pm  was 
a  possible  candidate.  The  current  work  suggests  that  1.6  or  2.3  pm  might  also  be  useful  choices, 
especially  if  several  such  channels  can  be  used  together.  Taken  together  several  VIS  and  NIR 
channels  can  be  used  to  provide  measures  of  both  cloud  geometry  (thickness)  and  microphysics 
(mode  radius  and  ELWC).  For  example,  a  shorter  wavelength  could  be  used  to  estimate  cloud 
thickness,  and  then  the  shortest  wavelength  unaffected  by  cloud  thickness  could  be  used  to 
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estimate  microphysics.  In  combination  with  other  data  or  other  satellite  sensor  information,  these 
reflectance  measurements  could  provide  even  more  accurate  estimates. 

Finally,  with  respect  to  ship  tracks  specifically,  the  results  above  should  have  some 
applicability.  The  tracks  are  quite  wide,  typically  a  few  to  twenty  kilometers,  compared  to  their 
depth,  typically  some  fraction  of  a  kilometer  (the  thickness  of  marine  stratus  and  stratocumulus 
clouds).  Thus,  the  assumption  of  an  infinite  plane  parallel  atmosphere  is  not  so  bad. 
Furthermore,  for  these  kinds  of  clouds,  changes  in  reflectance  are  probably  mostly  due  to  changes 
in  droplet  size  because  of  saturation  effects  in  the  other  possible  causes,  such  as  ELWC  and 
cloud  thickness.  The  model  results  show  that  reflectance  as  droplet  size  decreases  will  increase 
at  both  .63  pm  and  3.7  pm,  but  that  percentage-wise,  this  effect  is  more  pronounced  at  the  longer 
wavelength.  This  is  consistent  with  AVHRR  images  of  ship  tracks  in  channels  1  and  3.  Figure 
15  is  approximately  the  same  scene  as  in  Figure  1,  but  it  includes  radiance  scans  across  the 
tracks,  indicating  in  a  more  quantitative  way  how  the  radiance  is  enhanced  across  the  track.  It 
should  be  noted  that  channel  3  radiance  includes  both  solar  reflectance  and  thennal  emittance. 
When  this  thermal  component  is  removed  (by  estimating  the  Planck  temperature  from  channel 
4  (11  pm),  the  resulting  image  is  usually  very  similar  to  the  channel  3  image,  and  so  is  also  in 
agreement  with  the  model  results.  It  should  also  be  remembered  that  the  model  calculates 
hemispheric  reflectance,  whereas  the  AVHRR  images  depend  upon  sun-scene-satellite  geometiy. 
6.2  Future  Work 
6.2.1  N-layer  model 

Real  clouds  have  differing  droplet  number  concentrations  from  base  to  top.  If  the 
wavelengths  of  interest  are  short  enough  to  "see"  several  such  layers,  then  a  simple  one-layer 
model  may  not  be  accurate  enough.  Codes  for  simple  n-layer  models  are  readily  available  and 
should  be  tried  for  problems  such  as  cloud  models  with  an  evaporation  layer  on  top. 
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Figure  15.  Channel  3  scans  across  ship  tracks  in  the  eastern  North  Atlantic.  DTG  is  890512 
14:40  GMT.  Tire  scene  is  approximately  the  same  as  for  Figure  1. 
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6.2.2  Multi-channel  correlations 


Even  for  the  simple  model  used  for  this  report  much  of  the  data  base  remains  unexamined 
because  of  its  size.  Because  the  data  is  in  tabular  form,  it  should  be  straightforward  to  use  a 
database  management  system  to  look  for  other  variable  dependencies,  especially  involving  several 
channels  simultaneously.  This  offers  the  advantage  of  determining  some  parameter  in  several 
different  ways,  thus  improving  its  accuracy. 

6.2.3  Nonhomogeneous  finite  clouds 

Real  clouds  are  generally  finite  in  lateral  extent,  introducing  the  problem  of  side  effects, 
with  horizontal  composition  inhomogeneities.  To  model  such  clouds  requires  much  more  compli¬ 
cated  programs.  To  the  extent  that  it  is  possible,  keeping  a  model  as  simple  as  possible  is 
desirable,  since  it  makes  cause  and  effect  analysis  much  easier. 

6.2.4  Polarization  effects 

In  this  study  unpolarized  radiation  was  assumed.  It  would  be  interesting  to  investigate 
if  any  parameters  have  a  polarization  sensitivity. 

6.2.5  Comparison  with  field  data  -  satellite,  aircraft 

The  ultimate  test  of  model-based  conclusions  is  comparison  with  field  data.  For  the  case 
of  cloud  reflectances,  this  means  in  situ  measurements  of  cloud  properties  like  thickness,  droplet 
number  densities,  and  ELWC,  and  concurrent  satellite  measured  reflectances.  Some  data  is 
available  from  past  field  experiments  such  as  FIRE  and  relevant  planned  experiments,  such  as 
a  compoirent  of  ASTEX,  should  provide  additional  data. 
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